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Ultrafast Switching of Antiferromagnets via Spin-transfer Torque 
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Picosecond switching of the staggered antiferromagnetic order is shown to be realizable through 
spin-transfer torques from a short current pulse. The coupled dynamics of sublattice magnetization 
is mapped onto a classical pendulum subject to gravity and a driving pulse, where switching occurs 
if the pendulum acquires sufficient kinetic energy during the pulse to overcome the maximum of 
the effective gravity potential. The optimal switching scheme is explored through the dependence 
of switch angle and magnetic loss on the duration and strength of the current pulse. The physics 
discussed here provides a general route towards multi-functional THz applications via the spin- 
transfer torque in antiferromagnetic materials. 
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Introduction .—Ultrafast manipulation of magnetic 
states is a fundamental and lasting issue of magnetic de¬ 
vice design. Recently, attention has been increased in ul¬ 
trafast switching of antiferromagnets (AFs) using a laser 
pulse, in which the magnetic moments are first excited by 
a short pulse that lasts hundreds of femtoseconds, then 
magnetic switching takes place through inertia after the 
pulse is turned off [1-3]. This switching mechanism ex¬ 
ploits the giant exchange coupling between neighboring 
magnetic moments commonly found in AFs. Since the 
switching occurs on the picosecond time scale, which is 
orders of magnitude faster than the conventional switch¬ 
ing of ferromagnets, it opens a new avenue for ultrafast 
recording and processing of magnetically stored informa¬ 
tion using antiferromagnetic materials. 

However, to build a viable device using AFs, it is 
desirable to ask if the laser impetus can be replaced 
by an electrical current. A promising candidate is 
the spin-torque magnetoresistive random-access mem¬ 
ory (MRAM), where the magnetization is driven by 
the spin-transfer torque (STT) [4]. While spin-torque 
MRAM has only been realized using ferromagnets, recent 
progress in both experiments [5, 6] and theories [7-15] 
has pointed out the possibility of current-induced exci¬ 
tations of AFs. Based on a microscopic calculation of 
the electron scattering across a normal metal/AF inter¬ 
face [15], it has been demonstrated that a precessing stag¬ 
gered field pumps spin current into the adjacent normal 
metal, and vice versa : a spin accumulation impinging 
on an AF drives the coherent dynamics of the staggered 
field. The STTs discovered in AFs not only shed light on 
the mutual dependence between electron transport and 
magnetization dynamics, but also opens up the exciting 
possibility of electric control of AF devices. 

In this paper, we show that the staggered field of an 
AF can be switched within several picoseconds by a short 
current pulse with spin polarization perpendicular to the 
easy-plane. The switching process is first investigated by 
solving the coupled Landau-Lifshitz-Gilbert (LLG) equa¬ 


tions, and then further elucidated by an effective pendu¬ 
lum model. The material estimations are based on the 
widely studied room temperature AF insulator NiO, and 
the optimal switching is discussed by considering how 
the switch angle and magnetic loss depend on the pulse 
duration and STT magnitude. In addition, a THz nano¬ 
oscillator based on the perpendicular geometry is studied. 
The physics discussed here provides a general route to¬ 
wards multi-functional THz applications via the STT in 
antiferromagnetic materials. 

Dynamics .— For easy-plane AFs such as NiO and 
MnO, suppose the hard axis is z and the in-plane easy- 
axis is x. Scaling everything with frequency, we ex¬ 
press the out-of-plane anisotropy by cja < 0, the in-plane 
anisotropy by cu a > 0, and the Heisenberg exchange in¬ 
teraction by uje > 0. Within the macro-spin approxi¬ 
mation, the AF is characterized by two classical vectors 
mi and m 2 representing the magnetization of the two 
sublattices. Their dynamics is captured by the coupled 
LLG equations 

rh 1 = ujeTUi x m2 + uj a mi x x x mi 

+ WAmizZ x Trii + ctm\ x mi, (la) 
rh 2 = ujetu 2 x mi + oo a m 2x x x m2 

+ UAm 2z z x rn 2 + am 2 x m 2 , (lb) 

where a is the phenomenological Gilbert damping con¬ 
stant. We define the staggered field as t = (mi — m 2 )/2 , 
and the magnetization as m = (mi + m 2 )/2. By defini¬ 
tion, they are subject to the constraints: m • I = 0 and 
m 2 + £ 2 = 1. In the exchange limit, |m| \£\ so that 

£ 2 ~ 1, and thus £■£ ~ 0. The recombination of Eqs. (la) 
and (lb) gives the dynamics of m and £ 

rh = uj a x x (m x m + £ x £) + ujaz x Km + £ z £) 

+ ol£ x £, (2a) 

£ = 2c oe£ x m + uj a x x ( m x £ + £ x m) 

+ uo az x (m z £ + £ z m), (2b) 
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FIG. 1: (Color online) For p || z, an AF precession is im¬ 
plemented by the exchange torque, whereas a ferromagnetic 
precession resorts to the demagnetization field. 


where higher order damping terms like an x rh, am x h, 
and am x rh have been neglected since \m\ <C \£\- In 
Ref. [15], we derived the STTs that exert on m and £ in 
dimensions of frequency as 



y 


(a) 









(b) 




m m 


a 3 

Tn — z^G r £ X (£x V s ), 

(3a) 

a 3 

n = x ( m x V 8), 

(3b) 


where G r represents the real part of the spin-mixing con¬ 
ductance, V is the system volume, and V s is the spin 
voltage that impinges on the normal metal/AF interface. 
While Eqs. (3a) and (3b) are supposed to be added to 
Eqs. (2a) and (2b), they can equally well be decomposed 
into T\ — uj s mi x (pxmi) and r 2 = co s m 2 x ( pxm 2 ) and 
added to Eqs. (la) and (lb), where p is the unit vector of 
the spin polarization, and uj s = a 3 G r V s /(eV ) scales lin¬ 
early with the current density and inversely with the film 
thickness. When p || z, the two torques drag mi and m 2 
slightly out-of-plane as in Fig. 1, so that the exchange in¬ 
teraction generates precessional torque on the magnetic 
moments. It is this torque that switches the staggered 
field. Restricted by symmetry, the magnetization m de¬ 
velops only an out-of-plane component m = m z z. Corre¬ 
spondingly, the staggered field has only in-plane compo¬ 
nents so that i z — 0. These are confirmed by a straight¬ 
forward numerical simulation of Eqs. la and lb (see Sup¬ 
plementary Materials [16]). We mention in passing that 
no appreciable difference is observed in the switching be¬ 
havior between a compensated interface where both sub¬ 
lattices are subject to STT and an uncompensated inter¬ 
face where only one of the two sublattice is affected by 
the STT [16]. 

The system is now characterized by three variables 
{m z ^£ x ^y)- To formulate an effective description of the 
switching, we eliminate m z in terms of £. Taking the 
cross product of £ on Eq. (2b) leads to 


m = 


1 

2 uj e + \ua \ +u a £l 


£ x £. 


( 4 ) 


In typical transition metal oxides such as NiO, uo a and 
\uja\ are orders of magnitude smaller than uj e , thus we 
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FIG. 2: (Color online) Effective model of the staggered field 
switching, (a) perpendicular STTs cant mi ,2 slightly out 
of the basal plane, (b) by 0 — 20, the switching is mapped 
onto a pendulum subjected to gravity, damping, and a driving 
torque, (c) and (d): phase portraits of the pendulum for 
T = 0 and T > 0 (large enough to open up the attractor at 
origin), respectively. 


can safely ignore the last two terms in the denominator 
of Eq. (4). Substituting Eq. (4) into Eq. (2a) results in 


£x 


£ — u 2 r £ x x + 2auo E £ + 2uj e uj s (p x t) 


= 0, 


( 5 ) 


where ujr = yj2uj a uj E defines the frequency of in-plane 
oscillation in an AF [17]. Equation (5) describes the ef¬ 
fective dynamics of the staggered field £ in the exchange 
limit. If we focus only on the in-plane rotation of £, the 
degrees of freedom are reduced to one and Eq. (5) can be 
further simplified. Denote the in-plane polar angle by 0 
so that (£ x ,£ y ) = (cos 0, sin 0); then Eq. (5) gives 

uj 2 

4> + sin(20) + 2 auo E (j) = 2c d E u 8 p(t), (6) 

which is a damped non-linear oscillator with a driving 
force. Under the variable change 0 = 0/ 2, Eq. (6) is 
recast into 

0 + uj 2 r sin(6 > ) + 2auj E 0 = Aw e uj s p{t), (7) 

which is reminiscent of a pendulum subject to gravity and 
driving force as depicted in Fig. 2 (b). The pendulum is 
connected to a fixed pulley through a rigid but massless 
rod. If we twirl the pulley quickly by a short pulse, the 
motion of the system is described by JO + mgL sin 6 + 
aO = T(t), where J is the moment of inertia, mg is 
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the gravity of the ball, L is the length of the rod, a is 
the damping constant, and T(t) is the driving torque 
exerting on the pendulum. Upon the analogies: uje —> 
1/(2J), c o a mgL , and w s p(t) T(t)/ 2, Eq. (7) is 

mapped exactly onto the pendulum motion, by which the 
staggered field switching ( 7 r rotation of (j>) is represented 
by the crossing of the gravity maximum by the pendulum 
(27r rotation of 6). The phase portrait of the pendulum 
with (without) the driving torque T is plotted in Fig. 2 
2(d) [Fig. 2(c)]. When T = 0, the phase point (0,0) is 
an attractor; when T is sufficiently large, the attractor 
is broken up and flows towards the right, which enables 
the switching. 

During the pulse, energy is transferred from conduc¬ 
tion electrons to the sublattice magnetization m\ and 
m 2 via STT, and then stored in the exchange energy as 
they are canted non-collinearly. It is the releasing of this 
stored energy that provides an effective inertia to the mo¬ 
tion afterwards, through which the system surmounts the 
anisotropy barrier and finally relaxes to a new configura¬ 
tion. On the other hand, as shown in Fig. 1, a ferromag¬ 
netic switching under perpendicular spin polarization is 
engendered by the demagnetization field, which is orders 
of magnitude slower. While the exchange mechanism is 
intrinsic to AFs, the demagnetization field is sensitive to 
the shape anisotropy of ferromagnetic films. 

Switching .—With the effective pendulum model, we 
are able to perform a quantitative analysis of the stag¬ 
gered field switching. Using material parameters from 
NiO [17, 18], assuming a = 0.005 [19], we plot in Fig. 3 
the time evolutions of the out-of-plane magnetization 
2 m z ( t ) and the x-component of the staggered field l x ( t ). 
The switching process is characterized by the second of 
these, which is composed of two steps as follows: (1) the 
pulse drives the two magnetic moments slightly out-of- 
plane, thus they rotate under the exchange torque. (2) 
the staggered order moves to the opposite direction due 
to the inertia accumulated in the first step. The total en¬ 
ergy pumped into the system via STT finally dissipates 
away through the Gilbert damping. 

In Fig. 3, the canted magnetization during the switch¬ 
ing process is less than 1%, such that the approximation 
£ 2 « 1 is well respected. For small damping a = 0.002 
(not plotted), the flip is followed by a ringing tail; for 
a = 0.005 plotted in Fig. 3, the ringing effect is sup¬ 
pressed and the pattern resembles critical damping. Note 
that for small oscillations around the easy-axis, the crit¬ 
ical damping for NiO is roughly a crit ‘ = 0.0085; but for 
large angle rotations here, we observe a critical behavior 
for a smaller than a crit ‘. 

NiO has rocksalt structure, and all magnetic moments 
are attributed to Ni atoms that exhibit fee configura¬ 
tion. The {lll}-planes are ferromagnetically ordered and 
stack in an alternative manner, forming a layered AF. 
The easy spin direction is (112). To prepare the accu¬ 
mulation of conduction electron spins, we may adopt the 
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FIG. 3: (Color online) Out-of-plane magnetization (2 m z ) and 
staggered field projected on the easy-axis as functions of time 
in picoseconds. Pulse duration T — lOps; STT strength uj s — 
0.0034; Gilbert damping a = 0.005. Parameters for NiO: 
lje/( 27t) = 27.4THz, cj a /(27r) — 1GHz. 

spin Hall effect in a heavy metal [20, 21]. This method 
creates an in-plane spin polarization, thus the NiO needs 
to be grown in the {112}-direction to satisfy the required 
geometry in Fig. 1 and Fig. 2 (a). A much more effective 
spin accumulation can be achieved by driving the surface 
states of a topological insulator [22], which will signifi¬ 
cantly reduce the required current density. But in the 
following, numerics are restricted to the spin-Hall-driven 
NiO switching on a Pt/NiO interface. 

Attention should be paid to the proper choice of the 
STT strength. The staggered field will retrieve to the 
initial position if the STT is insufficient, and will over¬ 
shoot if it is too strong. For a Pt(20nm)/NiO(3nm) het¬ 
erostructure assuming no disorder and roughness on the 
interface [15], the uj s in Fig. 3 is converted to a current 
density of 6 ~ 7 x 10 7 A/cm 2 . The required STT for 
a proper switching becomes smaller when the pulse du¬ 
ration becomes longer. But for very long pulses, the re¬ 
quired STT ceases to reduce further. Therefore, to better 
understand the pulse dependence of the switching, we re¬ 
sort to the terminal angle of the staggered field and the 
total magnetic loss due to Gilbert damping as functions 
of the pulse duration and the strength of STT. 

As shown by Fig. 4, a smaller damping gives rise to 
a narrower window of the desired rotation (lowest or- 
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FIG. 4: (Color online) Upper left (right): terminal angle of 
l as a function of the pulse duration from 2ps to 20ps, and 
(jj s from 0.003 to 0.005 for a — 0.005 (a = 0.01). Lower 
left (right): total magnetic loss due to Gilbert damping for 
a — 0.005 (a — 0.01). 

ange region marked by i r). In real experiments, ampli¬ 
tude fluctuation of an electric pulse is inevitable, thus to 
stabilize the functioning of the device, extremely small 
damping is not favorable since it may easily lead to over¬ 
shoot. The total magnetic loss has a similar pattern as 
the terminal angle, hence it can hardly provide a prac¬ 
tically independent criterion for the optimal STT. How¬ 
ever, if Joule heating in the normal metal is taken into 
account, which scales as cj 2 T (T is the pulse duration), 
we obtain a different pattern compared to Fig. 4. To min¬ 
imize Joule heating with given magnetic loss, a shorter 
pulse with relatively stronger current density is preferred. 
Hopefully, with progress in nano-technology, picosecond 
or shorter current pulses can be realized in the future. In 
that case, the staggered field acquires a sufficiently large 
kinetic energy (angular velocity) long before it reaches 
the potential maximum of the in-plane anisotropy. The 
process would then be similar to the laser-pulse-induced 
AF switching [1-3]. However, a laser pulse couples to the 
AF via the weak Zeeman interaction, while STT rests on 
the strong exchange coupling between conduction elec¬ 
trons and local moments. 

Terahertz Oscillator. —Another possible applica¬ 
tion associated with the perpendicular geometry illus¬ 
trated by Fig. 1 is the terahertz nano-oscillator. If we 
replace the pulse current by a dc current in the Pt, the 
NiO will undergo a continual rotation. In Fig. 2 (b), it 
amounts to a continual rotation of the pendulum against 
gravity under a constant driving torque. Since uj a <C \ua\ 
in NiO, the in-plane anisotropy barrier causes only tiny 
non-uniformity in the angular velocity, which we can 
simply ignore. A continual rotation occurs when STT 


exactly compensates Gilbert damping. By a simple ge¬ 
ometry, the compensation condition is uj s = 2qluje sin$, 
where d is the azimuthal angle of i with respect to the 
hard £-axis. To leading order, the achieved oscillation 
frequency is linear in the STT (or the applied current): 

^rot. = bJ s /oi i ( 8 ) 

which is independent of the exchange interaction. For 
a = 0.002 and a Pt(20nm)/NiO(3nm) structure dis¬ 
cussed above, a ITHz oscillation requires a current den¬ 
sity of roughly 2^3x 10 8 A/cm 2 ; a thinner NiO and/or 
heavy metal with larger spin Hall angle will scale down 
the required current density. Moreover, the linearity of 
Eq. (8) remains robust up to about lOTHz, bearing an 
extraordinarily large range of linear control on the nano¬ 
oscillator. 

Just as with the nano-oscillator predicted by the au¬ 
thors of Ref. [15], the oscillation direction is determined 
by the current direction (a binary selection rule), which 
is unattainable in ferromagnetic materials. However, 
the excited modes are fundamentally different here. In 
Ref. [15], the spin accumulation in the normal metal is 
parallel to the staggered field, which excites the AF reso¬ 
nance eigenmodes. In this paper, the spin accumulation 
is perpendicular to the staggered field, and the switching 
is realized by exciting the spin superfluid mode [14]. 
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